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Abstract 

A three-dimensional anisotropic quantum well placed in an adiabatically precessing uniform mag- 
netic field is considered and an explicit formula for the Berry phase is obtained. To get the Berry 
phase, a purely algebraic algorithm of reducing a quadratic Hamiltonian to the canonical form via 
symplectic transformations of the phase space is presented. 
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I. INTRODUCTION 

n 

In the classical work M. Berry has shown that if an evolution of a physical system is 
cyclic and adiabatic the geometric phase (so-called Berry phase) is accumulated on the wave 
function in addition to the usual dynamical phase. The Berry phase has a purely geometric 
nature and reflects geometric properties of the parameter space of the system. B.Simon 
has given an interpretation of the Berry phase as due to holonomy in a line bundle over the 
parameter space. 

The space of values of the magnetic induction as a parameter space for the considered system 
is frequently dealt with. As a rule, for a magnetic field, the Berry phase is studied in the two 
cases: (1) a potential well is transported along a closed loop in the usual configuration space 
R U, M (2) spin particles interact with a time-depended magnetic field ((J The most 

part of theoretical studies is devoted to the evolution of the spin in a magnetic field; in this 
case explicit formulae are easy to obtain using finite-dimensional algebra. However, the Berry 
ohase can arise in the case of a spin-less particle interacting with a precessing magnetic field 
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The aim of this paper is to obtain an explicit formula for the Berry phase in the case of a 
three-dimensional anisotropic quantum well placed in a precessing magnetic field. The classical 
Hamiltonian of such a system is a quadratic form of momenta and positions. The considered 
problem has attracted attention of physicists since the 1980's. In particular, M. Berry jl^] 
has calculated the geometric phase for a one- dimensional quadratic potential with variable 
coefficients. As for now, the most general result for quadratic Hamiltonians is contained in j^| 
where the following simple case has been studied: a quadratic potential has an axial symmetry 
with respect to the z-axis and a magnetic field precesses around z-axis. These assumptions 
reduce the dimension of the problem and significantly simplify calculations. We investigate 
the most general problem in the case in which the frequencies of the parabolic potential are 
mutually different and the magnetic field is arbitrarily directed. 

The significant peculiarity of a quadratic quantum Hamiltonian lies in the fact that the 
Feynman integrals for their propagators can be calculated explicitly. What is more, the problem 
of calculating these integrals is reduced to finding the canonical form for a classical counterpart 
of the Hamiltonian via symplectic transformations of the phase space; in our case this form is 
a sum of squares of momenta and positions). Finding the canonical form is a standard (but 
as a rule very cumbersome) problem of linear algebra. For convenience of readers, in the next 
section we present an algorithm of finding the classical (symplectic) and quantum (unitary) 
transformations that reduce the corresponding Hamiltonian to the canonical form. 
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It is important to stress that we obtain our results using only simple methods of linear 
algebra. It is strange enough that in the similar situations, more complicated methods are used 
in the physical literature; in particular, the method of the Bogolubov canonical transformations 
or other indirect methods. Note another adva„ta g e of the method we use. To find various 
physical quantities one needs the matrix elements of the corresponding operators. However, 
a direct calculation of matrix elements can be a complicated problem. Our method allows to 
resolve this problem using only simple calculations from linear algebra. As an example we 
can mention the papers [12J and where the hybrid and hybrid-phonon resonances for an 
asymmetric quantum well in a magnetic field have been studied by means of this approach. 



II. PRELIMINARIES 



The advantage of a description of mechanical systems with the help of the Hamilton for- 
malism is conditioned by the following main reason: a choice of the generalized coordinates is 
restricted by no conditions (any one of quantities defining the state of a system may be the 
generalized coordinates.) What is more it is possible to choice such a class of transformations 
of phase coordinates that includes the transformation of the all initial phase coordinates pi, qi 
to the new ones Pi, Qi (here i — 1, 2...n; n is the dimension of the configuration space) because 
momenta and positions play the same role in Hamilton equations 

Q = Q(p,q);P = P(p,q). (1) 

However there is a restriction. The transformation is permissible (canonical) if and only 
if it conserves the Poisson brackets. The equivalent criterion for the transformation to be 
canonical is the equality to the exact differential the 1-form pdq — PdQ. This 1-form is the 
exact differential of the generating function of the canonical transformation (£[]). 

dS = pdq - PdQ. (2) 

Note that if we know the generating function we can always find the functions determing the 
transformation (|T|). 

Through the paper we use the notation 

fciH = -(pi|fe) + (P2\qi) (3) 

for the standard symplectic product in the phase space M 2ri = R" x M2 ; here Xj = {pj,qj), 
(j = 1,2). 
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Recall that a linear transformation 

n \ / p \ 

(4) 

is canonical (i.e., preserves the Poisson brackets) if and only if the matrix L is symplectic: 

[Lxi\Lx 2 ] = [xi\x 2 \. (5) 
It is convenient to write the matrix L in the form 

(6) 

The following commutation relations between components of the matrix L follow immediately 
from the fact that the matrix L is symplectic 

L\L/± — L2L3 = E , L3L4 — L4L3 = , L\L 2 — L* 2 Li = , (7) 

where E is the unit matrix. 

Denote by S the generating function for the transformation L (this means by definition that 
dS = pdq — PdQ); we will suppose that S is a function of p and P. Then 

dSjjp, P) dS{p, P) f 

Q = — dp-' q = ^p~- (8) 

In this case the generating function of the canonical transformation Q is expressed in terms 
of Li by the following formula (we assume that I-L2I 7^ 0) 

S{p,P) = -(P\L?p) + i(P|L 2 %P> + ^(p\L 4 L^p). (9) 

According to the Dirac-Fock theorem [14| the classical transformation (J3| corresponds to a 
quantum transformation U : L 2 (Rp) — ► L 2 (]R™) that is to a unitary operator related to L as 
follows 

U- l pU = L X P + L 2 Q, 

(10) 

U- 1 qU = L 3 P + L 4 Q. 
The operator U can be represented as an integral operator 

(Uf)( P )= [ U(p,P)f(P)dP (11) 



with the kernel 



U(p,P) = Cexp 



~S(p,P) 



(12) 
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where the normalization constant C = ((2nh) 3 \ det Z^l)" 1 / 2 . Then using the unitary of the 
operator U we can express the initial wave function of the quantum Hamiltonian in terms of 
the final wave functions by the following formula 



ij;(p) = C / ^(P)exp 



- h S(p,P) 



dP, 



(13) 



here ip(p) is the wave function of the original Hamiltonian and ip(P) is the wave function of 
the unitary equivalent Hamiltonian. 



III. DIAGONALIZATION OF QUADRATIC HAMILTON! ANS 

The purpose of this section is to demonstrate an algorithm of reducing the quadratic Hamil- 
tonian 

-<»■>-*(:)*(:)• 

where if is a symmetric matrix of order 2n x 2n, to the canonical form. This form must contain 
only squares of the corresponding momenta and positions. 

Then Hamilton equations corresponding to our quadratic Hamiltonian (pPfj) can be written 
in the form x = IHx, here I is the symplectic unity 



—E 
E 



(15) 



At first we find the family (Aj)i<j<2 n of eigenvalues the matrix IH, where I is the symplectic 
unity 

(o-e\ 

I=\ . (16) 

and the family (fi)i<i<2 n °f the corresponding eigenvectors. For this purpose we use the fol- 
lowing simple lemma 

Lemma. The characteristic polynomial det (7 H — XE) has only real coefficients and contains 
only terms of even degree. 

It follows from the Lemma that if Aj is an eigenvalue of the matrix IH then — Aj and A* are 
eigenvalues of this matrix too. Taking into account this assertion we can arrange Aj in a finite 
sequence of the form 

u u .., v h iui, ..iu m , -v u .., -v h -iu u .., -iu m , (17) 
where z/j, uj. > 0; I, m > 0; I + m — n. 
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Now we transform the basis (/j) into a symplectic one. Let (&j)i<j<2n be the following 
sequence of vectors from R 2 ™: 

fi(vi),~, fi(yi), Re[/i(«Ji)],.., Re[f m (iu m )), fi(-vi), .., fi(-vi), lm.[f 1 (iuj 1 )), .., Im[/ ffl (iw m )]. 

(18) 

It is clear that the vectors (&i)i<j<2n form a basis of the space R 2n . These vectors satisfy the 
skew-orthogonality conditions 



[bi\bk] =0, if i < k; i + n ^ k, 
[bi\b i+n ] ^ 0, V i = 1, ..,n. 

Denote q = [6j|6 i+n ], then according to (fT9|l the vectors b[, 



(19) 



sgn Ci . 

6, = 7p=rOi; « = i, 



6' 



(20) 



i+n 



'■ bi+m i 1) ■ ■ ; 



form a symplectic basis of the space R 2n . Using coordinates with respect to this basis, we 
reduce the Hamiltonian H to the form 

l m 



2 ) 



(21) 



i=l i=l 

where Si = sgn q. 

Now we define the canonical basis of the quadratic form H by the relations 

di = b[, if I + 1 < z < n and n + I + 1 < z < 2n, 
fl i = -^( 6 i + 6 i+»). if !<»*<'. 



(22) 



With respect to this basis, the Hamiltonian H has the canonical form 

l m 

h(p, Q) = Y< v < ( p ? -Qt) + J2 £ ^ + Qli) ■ ( 23 ) 

i=l i=l 

The transition matrix from the initial phase coordinates p, q to the final ones P, Q can be 
written in the following form 

L — [ai, a 2 , a 3 , a 4 , a 5 , a 6 ] , (24) 
where each vector a, is written in the corresponding column. 
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IV. DIAGONALIZATION OF THE HAMILTONIAN OF A 3D ANISOTROPIC OS- 
CILLATOR WITH A MAGNETIC FIELD 



Let us consider the Hamiltonian of a 3D anisotropic harmonic oscillator with a magnetic 
field B arbitrarily directed with respect to the potential symmetry axes 



(25) 



Here A is the vector potential of the magnetic field, c is the light velocity, flj (i = x, y, z) are 
the characteristic frequencies of the parabolic potential, m, e, and p are the mass, the charge, 
and the momentum of the considered particle. 

The results of Sec. Ill show that Hamiltonian (J2*H|) is unitary equivalent to a quadratic 
Hamiltonian without magnetic field but with new frequencies named "hybrid frequencies". 
The purpose of this section is to reduce Hamiltonian (J23)) to the canonical form 



ff(P, Q) = ^ ( p i + Pi + Pi) + 7 HQl + "IQI + "IQD 



(26) 



using the algorithm developed in the preceding section; here P, Q are the new phase coordinates, 
LUi (i — 1, 2, 3) are the hybrid frequencies. It is convenient to choose the following gauge for the 
vector potential A 

A = 



^B y z - B z y, 0, B x y - ^B y x ) . 



(27) 



Below we use the notations 



K\ = m 2 Q 2 x , K\ = m 2 Q 2 K 2 = m 2 fl 2 z 



Px = Pi, Py = P2, Pz =P3', x = q u y = q 2 , z = q 3 ; 



eB x 
2c 



B X , ^ 
U 2c 



B 2 , 



2c 



B,. 



In this notation 



H = 2^ ~ B2Q3 + 2Bsq2)2 +P2 + (Ps + B2Ql ~ 2Bl<l2)2 + R2q2 + K2q2 + K2q ^ ' (28) 



i.e., H is a quadratic form of p and q with the matrix 



H 



1 

rn 



/ 


1 











2£? 3 


-B 2 


\ 







1 
























1 


B 2 


-2B X 















B-2 


Ki 


-2B X B 2 









2B 3 





-2B X 


-2B X B 2 


K 2 


—2B 2 B 3 




\ 


-B 2 











—2B 2 B 3 


K 3 


J 



(29) 
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where K x = K\ + B 2 , K 2 = K 2 + AB 2 + AB 2 , K 3 = K 2 + B\. 

To find the eigenvalues of the matrix IH we solve the equation det(IH — XE) = 0, which is 
easily reduced to the form 



1 
1 


n 
u 


n 
u 


-fl 


9 R 


— E>2 


n 
u 


i 


n 
u 


n 
u 




n 
u 


u 





1 


B 2 


—IB\ 


-fi 











K\ + f 


-2B 3f x 


2B 2f i 











2B 3f i 


K 2 + f 


-2B x n 











-2B 2f i 


2B xl i 


K\ + f 



(30) 



where /i = mX. Since 




det(AC) 



(31) 



for any nxn matrix A. B,C, we can obtain the hybrid frequencies from the following sixth-order 
algebraic equation [15] 

+ A 2 )(fi? + X 2 )(Q 2 Z + A 2 ) + ulM + A 2 )A 2 + u 2 yc (Q 2 y + A 2 )A 2 + u 2 zc (Q 2 z + A 2 )A 2 = 0, (32) 

where Ui C = eBi/mc are the components of the cyclotron frequency and A^ = ±iuji. Note that 
equation (J3*2j) always has three different real negative solutions with respect to A 2 . 

Now our purpose is to find the eigenvectors of the matrix IH associated with the eigenvalues 
A^. Using elementary matrix transformations one can reduce the matrix M = IH — XE to the 
form 

-K\-ii 2 2B 3f i -B 2l i \ 
H 2 2B xl i 
-2B xi i 
2B 3 

-2B X 



m 



/o 





1 
1 

\0 1 



-2B 3l i 
2B 2f i 
-// 


B 2 



-Kl 



-B 2 





(33) 



The eigenvectors of ()33|) are obtained from the system 



(34) 



Xq = 0, 
P = -^q- 

Solving this system we get the coordinates of the eigenvectors g(fi) associated with the eigen- 
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values [i = m\ i of IH: 

Pl = -2B 3 K\{K 2 + /2 2 ) - 4B 2 B 3l 2 2 - 2B 1 B 2l 2(K 2 - /i 2 ), 
p 2 = ^K 2 + ^ 2 )(K 2 + ^ 2 )+AB 2 ^, 
p 3 = 2B 1 K 2 (K 2 + fi 2 ) + 45 1 5 2 V - 2B 2 B 3f i(K 2 - fi 2 ). gi = 2B 3f i(K 2 + /i 2 ) + 4B 1 B 2f i 2 , 

q2 = {K 2 + ^ 2 ){K 2 + ^ 2 )+AB 2 ^ 2 , 
q 3 = AB 2 B 3f i 2 - 2 f iB 1 (K 2 + fi 2 ) . 

(35) 

Consider the vectors 

fi = Re \g(-imwi)], f 2 = Re \g(-imw 2 )\, f 3 = Re [g(-imu> 3 )], 

/ 4 = Im [g(-imwi)\, f 5 = Im [g(-imw 2 )], f 6 = Im [g(-imw 3 )\. 
Then we have for the coordinates of fi the following expressions. If i — 1, 2, 3, then 

Repi = -2B 3 K 2 (K 2 - m 2 cu 2 ) + 4 J B 2 J B 3 m 2 cc> 2 , Rep 2 = 0, 

Rep 3 = 2B X K\{K\ - m 2 u 2 ) - AB^im 2 ^ 2 , 
Regi = -45i B 2 m 2 uj 2 , Req 2 = [K\ - m 2 ^ 2 )^ 2 - m 2 cu 2 ) - 4 J B 2 m 2 w 2 , Reg 3 = -AB 2 B 3 m 2 uj 2 . 
If i = 4, 5, 6, then 

Inrpi = 2B 1 B 2 muj i ^ 3 (K 2 +m 2 uj 2 __ 3 ), lmp 2 = ~muJi^ 3 (K 2 -m 2 uj 2 _ 3 )(Kl-m 2 cj 2 _ 3 ) +AB 2 m 3 ujf_ 3 , 

lmp 3 = ^B^mUi^K 2 + m 2 u 2 _ 3 ), 
lmq 1 = -2B 3 muJi_ 3 (Kl - m 2 uj 2 _ 3 ), lmq 2 = 0, Img 3 = 2muj i _ 3 B 1 (K 2 - m 2 oo 2 _ 3 ). 
Denote 

M t = m* {co^l + u 2 yc ) + uljtfM - + - c 2 ) 2 + 



+ [(^-^ 2 )(^-^ 2 )-<c 2 ] 2 } 1/2 (36) 

for i = 1,2, 3; and = M;_ 3 for % — 4, 5, 6. It is easy to show that the following basis is 
symplectic 

fH = 1 ft , (37) 

where we put uJi = cjj_ 3 for i = 4,5, 6. 

With respect to this basis the Hamiltonian H has the form 

3 

2 



# = *X>(# + #)> (38) 



i=l 
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here pi, and qi are new phase coordinates. 

To reduce this Hamiltonian to the canonical form (J2*fi|) we change the variables once again 

Pi = y/muipu Qi = , Qi- (39) 

Finally, the transition matrix L from the initial coordinates to the coordinates Pi, Qi has 
the components 

lij = M~Jbr (40) 

where defined as follows: 



if i — 1, 2, 3 then 

5 n 2/n2 2\ i 1 5 2 2 n 5 r>2/r>2 2\ ^ 5 2 2 

Oil = -m u^O^O^ - ) + -m u)^^ , a i2 = 0, a i3 = m uj xc il z [il x -u i )--m u) X( w yc u) i 



a iA = -m A uj xc uj yc uj 2 , a i5 = m 4 (VL 2 x - u 2 )(Q 2 z - uj 2 ) - m 4 u yc uj 2 , a i6 = -m 4 u yc u zc uf; 
if i — 4, 5, 6 then 

a a = ^m 5 u xc u} yc u}i_ 3 {Q 2 x + w?_ 3 ), a i2 = -m 5 ^_ 3 (fi 2 _ w ?_ 3 )(^ - w?_ 3 ) + m 5 u 2 yc ujf_ 3 , 

«i3 = ^m 5 uj yc uj zc uj i _ 3 (Q 2 z + uj 2 _ 3 ), a iA = -m 4 w 2C u; i _3(fi^ - uj 2 _ 3 ), 

a i5 = 0, a i6 = m A u xc Ui-z(VL 2 x - uj 2 _ 3 ). 

As a result we have obtained an explicit formula for the transition matrix L. 

Note that the quantum counterpart of classical Hamiltonian (121)1 is the quantum harmonic 
oscillator with the eigenvalues 

E nin2n3 = fkoi \ ni + ~~ \ +tw 2 (ni + 0+^3 (n\ + ^ , (41) 
where ri\, n 2 , n 3 = 0, 1, ... In the P-representation, the corresponding eigenfunctions are 

MP) = {U^y 1 ^ X n 2 Xn 3 (~p±) • (42) 

Vv^ift/ Vv w iV Vvwift/ 

Here Xn is the n-th oscillator function 

X n(x) = (vr 1 / 2 2"n!)- 1 / 2 exp(-x 2 /2)P n (x), 



where H n is the corresponding Hermite polynomial. 
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V. CALCULATING THE BERRY PHASE 



The purpose of this section is to obtain an explicit formula for the Berry phase of a 3D 
anisotropic parabolic quantum well placed in a precessing magnetic field. We deal with the 
case when a considered quantum state undergoes an adiabatic evolution in the space of the 
values of the magnetic induction B in such a way that eigenvalues (|41j) remain nondegenerate. 
Let ip n be the wave function of the original state (here n = (ni, n 2 , n 3 ) is the set of the quantum 
numbers). If the system undergoes an evolution in question, then the wave function obtains 
the geometric phase 

7n = y"v n (B)dB, (43) 

G 

where C is a closed path B(i) in the parameter space such that B(T) = B(0), and V n (B) is 
so-called Berry vector potential given by 



V n (B) = i(^(B)|V B Vn(B)) = -Im(^(B)|V B ^(B)>. 



(44) 



Since the Berry phase is determined by V n (B), it is sufficient to calculate only the vector 
potential. In our case the magnetic field doesn't vanish and hence the parameter space is 
R 3 \{0}, which is topologically non-trivial, and we can expect that the Berry phase is non- 
trivial, too. 

To find the scalar product in (pf4*|) one needs eigenfunctions of the initial Hamiltonian. How- 
ever, a direct calculation of the wave functions to obtain the Berry phase is a complicated 
computational problem. To calculate the Berry phase, we suggest to simplify the calculations 
using the method of linear canonical transformation of the phase space developed in the preced- 
ing sections. In this case we can use formula (|TT|) to express the initial wave functions ^(P; B) 
in terms of the final ones 



^ n (p;B) = C / ^ n (P;B)exp 



h 



S(p,P) 



dP. 



Taking into account that ^(P; B) is real-valued we can calculate VeV ; n(p; B): 



V B ^ n (p;B) 



V B Cty n (P;B) 



cxp 



-^(p,B) 



dP- 



(45) 



-j: J C^ n (P;B)exp 

In what follows we denote V B simply by V. 
Let us calculate scalar product 



-^(P,P) 



(^n|V^n)) = / d P <(p)Wn(p) 



[V B 5( P ,P)]dP. 




dp<C* / dPdPVn(P') VC^ n (P 



(46) 
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x exp 



h 



(5(p,P')-5(p,P)) 



h 




dPdP'|C|Vn(P>n(P)x 



x exp 



-(5(p,P')-S(p,P)) 



v^(p,p; 



(47) 



By changing the variables L 2 l p/h = y; dp — | det L 2 \h 3 dy, we perform the integration with 
respect to p in the first term of (|4*7)l and get 



exp 



h 



(5(p,P , )-5(p,P)) 



dp 



exp 



^((P'\L 2 %P'> - (Pi^-^xP)) 



det L 2 \h 3 J exp [i(P - P'|y)] rf y 



= (27rft) 3 | det L 2 \5(P - P'). (48) 

Hence the first term in ()47|) is real- valued and doesn't contribute in the Berry phase. Therefore, 
we need to calculate only the second term in (|47|). 
It is convenient to introduce the notations 

A = -VL^ 1 ; B = l -V(L 2 %); D = ^(L^J 1 ); A = ^((P'| J LJ 1 J L 1 P')-(P| J LJ 1 L 1 P)). (49) 

Now integrating the second term in ()47|) with respect to p, we obtain 



exp 



ft 



(5(p,P , )-5(p,P)) 



VS(p,P)dp 



ft 3 1 det L 2 1 exp (iA) x 
3 

(2tt) 3 (P| j BP5(P - P')> - ^(2tt) 3 J2i(AL 2 )*I>} k d k 6(P - P')- 



fc=i 



ft 2 (27r) 3 ((d|L*I}L 2 d))5(P - P')] ■ 
Calculating the corresponding derivatives 



<9 fc exp (iA) = -- [(L2 X LiP)J exp (iA) 



(50) 



(51) 



d^exp (iA) = -^(L2%)«exp (iA) - ^(L 2 ^iP)^ %P), exp (iA) , (52) 
ft ft 

and substituting (JHTJ) and ijH^jl into (JHUj) after simple algebra we get, taking away the imaginary 
part, 

ReJ = 
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(2nh) 3 \ detL 2 \ {{P'\BP'} + (P'lA^P') + (P'^DLiP') - h 2 (d\L*DL 2 d}} 8(P - P'). (53) 
As a result we have the following formula for the Berry vector potential 

V n (B) = -Im / dpC(p)V^ n (p) = \\ \ rfPrfP|C| 2 ^n(P , )^n(P)(Re/) = 



h 



(2nh) 3 \detL 2 \\C\ 2 y dP'[(P'\BP') + (P'\(AL 1 )P') + (P^DLOP)] |^(P')| 2 + 

3 

ft 2 / tZPcZP' ]T {LlDL^duSiP - P')^ n (P')^n(P) 



it follows for 7^ Z 



If k = Z, then 



Finally we have 



V n (B) = J2 (»* + 



(55) 



Denote 

F = .B + AL X + L\DL X , 
G = L* 2 DL 2 . 

It is clear that G* = G and we can rewrite (|54jl in the form 

V n (B) = ±(2nh) 3 \detL 2 \\C\ 2 {J dP / (P / |FP / )|^ n (P , )| 2 dP , + 

h 2 J dPdP' £ G^dkSp - P , )^„(P')^n(P) | • (56) 

Using (|4"2|) . we get 



ii k ^ I, and 



Since 



P fc P^„(P)| 2 dP = 0, 

y p 2 ^ n (p)i 2 tzp = ^ y ximm 2 ^ = + 1) ^ . 

Xn(x)Xni x ) dx = °> ( 57 ) 



V„(P)M^n(P)dP = 0. (58) 

y ^ n (p)9 2 ^„(p)rfp = -JL f nfc + . (59) 



Wfci 5 ** + — I . (60) 



Equation (|50|l is the main result of our paper. This formula extends a result of Berry to 
the 3D case and contains as a particular case some results from P|. 
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Remark. Note that Formula (|6U|) is valid in more general situations. Namely, let H(£) be a 
family of quadratic Hamiltonians in the state space L 2 (M 3 ) depending on a parameter (, ( 6 I, 
where the parameter space X is a smooth manifold. Let £ = be a curve in X such that 
H(£(t)) has three different eigenvalues at each t. Then the Berry phase corresponding to an 
adiabatic evolution along the curve £ = £(t) is given by (JHDj), if we replace Vb by in all the 
auxiliary expressions. 

Acknowledgments 

The authors are very grateful to J.Briining and K.Pankrashkin for useful discussions. We 
gratefully acknowledge grants of DFG, INTAS and RFBR (Projects No. 01-02-16564 and 02- 
01-00804). We are also thankful to Humboldt University of Berlin for warm hospitality during 
the preparation of this paper. 



[1] M. V. Berry: Proc. R. Soc. London, Ser A 392 (1984), 45. 

[2] B. Simon: Phys. Rev. Lett. 51 (1983), 2167. 

[3] M. V. Berry: J. Phys. A: Math. Gen 30 (1997), 8355. 

[4] S. -L. Zhu: Solid State Comm. 30 (2000), 233. 

[5] P. Exner, V.A.Geyler: Phys. Lett. A 276 (2000), 16. 

[6] T. Bitter, D.Dubbers: Phys. Rev. Lett. 59 (1987), 251. 

[7] D. J. Richardson, A. I. Kilvington, K. Green, S. K. Lamoreaux: Phys. Rev. Lett. 61 (1988), 2030. 

[8] H. Weinfurter, G. Badurek: Phys. Rev. Lett. 64 (1990), 1318. 

[9] D. J. Fernandez, M.A. del Olmo, and M.Santander: J. Phys. A: Math. Gen 25 (1992), 6409. 

[10] M. V. Berry: J. Phys. A: Math. Gen 18 (1985), 15. 

[11] H. A. Fertig, B. I. Halperin: Phys. Rev. B 36 (1987), 7969. 

[12] V. A. Geyler, V. A. Margulis, A. V. Shorokhov: Phys. Rev. B 63 (2001), 245316. 

[13] V. A. Margulis, A. V. Shorokhov: Phys. Rev. B 66 (2002), 165324. 

[14] P. A. M. Dirac: Lectures on quantum mechanics, Yeshiva University, New York 1964. 

[15] Q. P. Li, K. Karai, S. K. Yip, S. Das Sarma, H. D. Drew: Phys. Rev. B 43 (1991), 5151. 



14 



